kurtosis based Wavelet Filter
Chiarelli Antonio Maria
Beckman Institute of Technology
latest version 08/21/2014
Clears Motion artifacts from NIRS data by applying WK algorithm
Chiarelli, A. M., et al. (2015). "A kurtosis-based wavelet algorithm for motion artifact correction of fNIRS data." Neuroimage 112: 128-137.

inputs: signal: matrix, n time points x m channels, th_kurt: kurtosis
threshold (suggested 3.3), level: level of decomposition to start
applying the algorithm (>3=)
outputs: artifact cleaned signal: matrix, n time points x m channels